Fermi-liquid theory of imbalanced quark matter 
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The temperature dependence of the thermodynamic potential of quantum cliromodynamics 
(QCD), the specific heat, and the quark effective mass are calculated for imbalanced quark mat- 
ter in the limit of a large number of quark flavors (large-A^i?), which corresponds to the random 
phase approximation. Also a generalization of the relativistic Landau effective-mass relation in the 
imbalanced case is given, which is then applied to this thermodynamic potential. 
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I. INTRODUCTION 



Landau Fcrnii-liquid theory has seen considerable suc- 
cess in describing a wide variety of fermionic many- 
particle systems, such as liquid Helium-3, electrons in 
metals, nuclei and nuclear matter [l|-lj]. It gives an effec- 
tive description of the low-lying elementary excitations 
(quasiparticles) at low temperatures, i.e., at a tempera- 
ture T such that the system can be considered degener- 
ate (T <?; ^, where ^ denotes the chemical potential) but 
still in the normal phase above any symmetry breaking 
phase transition (T ^ Tc), for example, to a magnetic or 
superconducting phase. Nevertheless, the theory is not 
only important for the description of the above "normal" 
(Fermi-liquid) phase, but is also vital to correctly de- 
scribe the emergence of a possible ordered phase [5| . In- 
deed, according to the Bardccn-Cooper-Schrieffer (BCS) 
theory the onset of superconductivity is to be viewed as 
an instability of the Fermi liquid under an attractive in- 
teraction, which results in the formation of Cooper pairs 
[6| . The same is true for magnetism in the case of repul- 
sive interactions. 

In particular, BCS theory implies that the high-density 
and low-temperature region of the phase diagram of 
quantum chromodynamics (QCD) contains a color super- 
conducting phase of quarks [7|, |8|. In fact, since quarks 
carry color, flavor and spin quantum numbers, many dis- 
tinct superconducting phases arc possible and are char- 
acterized by the various symmetries of the Cooper-pair 
wavefunction. These phases of cold and dense QCD 
might occur in the core of neutron stars where mat- 
ter is compressed to several times the nuclear density 
po — 0.16 fm~ . The presence of such a superconducting 
phase is expected to have observable consequences on the 
cooling and magnetic fields of neutron stars [9l-[ll|. A 
particularly interesting phase is the so-called 2SC phase 
in which only two flavors of quarks are paired, while the 
unpaired quarks of the third flavor remain a Fermi liquid. 
This phase further illustrates the importance of under- 
standing the quark-gluon plasma, i.e., the normal phase 
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of quarks. 

Fortunately, due to asymptotic freedom, the coupling 
constant g of QCD becomes small at large chemical po- 
tential fi, such that a systematic study of the high-density 
and low-temperature region of the QCD phase diagram 
is possible using perturbation theory. At high densities 
the dominant interaction between quarks is that of one- 
gluon exchange, where the long-range behavior of the 
gluons is screened due to the quark-gluon plasma. While 
the electric gluons (longitudinal) are screened by a De- 
bye mass rrig ^ gfi, the ma gnet ic (transverse) gluons are 
only dynamically screened jl3 . [l3j . This residual long- 
range behavior of the magnetic gluons dominates the low- 
temperature behavior of the system. Examples of the ef- 
fects of magnetic gluons can be seen in the non-standard 
scaling of the critical temperature of the superconduct- 
ing phase Tc ^ fj.g~^ exp{-3TT^/g) [Tl - [l6J . which de- 
pends exponentially only on 1/g instead of showing the 
expected l/g^ behavior, or of the life-time of elementary 
excitations of (grand-canonical) energy E near the Fermi 
surface t{E) ~ \E\~^ |17H2g | that signal a (marginal) 
non-Fermi-liquid behavior at zero temperature. 

In the context of a possible deconfined phase of quarks 
inside neutron stars, the Fermi system is expected to be 
imbalanced due to the different chemical potentials of the 
various quarks. This deconfined phase may contain u, d, 
and s quarks with different densities as a consequence 
of their different masses and charges. In cold atomic 
physics, an analogous population imbalance has been re- 
alized experimentally between two spin states [2l|, |23 |. 
which resulted in frustration of the pairing between par- 
ticles. This will cause very different behavior from the 
normal BCS case and is currently a very hot topic in a 
wide variety of fields [2j, [2^. In particular, the imbal- 
ance between different flavors of quarks is expected to 
significantly alter the properties of the QCD phase di- 
agram, such as the 2SC phase mentioned above or the 
color-flavor locked phase where all the colors and flavors 
are paired [7|. In the present paper, we therefore gener- 
alize the quasiparticlc properties based on Fermi-liquid 
theory to an imbalanced system. 

To do so, the effect of the dressed gluons on the ther- 
modynamic potential, specific heat and effective mass 
of the quasiparticles is determined at low temperatures. 



The thermodynamic potential will be calculated for a 
two-flavor quark system in the limit of a large number 
of quark flavors (large- A^f)- This approximation cor- 
responds to the random phase approximation (RPA), 
which has been quite successful in condensed-matter 
systems such as the interacting electron gas in metals 
[25| . Using the framework of Landau Fermi-liquid the- 
ory for relativistic systems, derived in section |TT1 the ef- 
fective mass of the quasiparticles can be determined from 
the thermodynamic potential by considering the specific 
heat. The result and its implications on the applicabil- 
ity of Fermi-liquid theory is discussed. As we will show, 
the logarithmic dependence of the temperature, which 
is due to the transverse gluons, shows the breakdown of 
Fermi-liquid theory at low temperatures. Furthermore, 
even though within our approximations there is no inter- 
action between the different flavors of quarks in the limit 
of weak coupling, the RPA correction still presents some 
mixing between the different flavors, which stems from 
the long- wavelength screening of the longitudinal gluons. 
Above and throughout the article natural units are used, 
i.e., units such that /i = c = fcs = 1. Other conventions 
and technical details of the calculation can be found in 
the Appendices. 



II. LANDAU FERMI-LIQUID THEORY 

To describe the normal state of the quark-gluon plasma 
at nonzero temperatures, which is a strongly interacting 
gas of quarks, Landau Fermi-liquid theory can be used. 
This theory takes as a starting point the non-interacting 
Fermi gas and switches on the interaction adiabatically. 
As long as the temperature of the system is much higher 
than the critical temperature Tc for superconductivity, no 
bound states (Cooper pairs) will form and each state of 
the non-interacting Fermi gas is transformed into a state 
of the interacting gas. Therefore, the excitations of such a 
Fermi liquid remain of a fermionic nature. For simplicity, 
Fermi-liquid theory will be introduced here for a system 
at zero temperature even though the quark-gluon plasma 
is formally a marginal Fermi liquid in that case. However, 
the introduced concepts turn out to be valid for nonzero 
temperatures much lower than the Fermi temperature 
Tp, which is expected to be around 10^'^ K for quarks in 
the core of neutron stars. 

The basis of Landau Fermi-liquid theory is to consider 
the effect of a small change of the ground-state Fermi dis- 
tribution on the thermodynamic potential density, which 
can be written as 



5^ 



^n<,d><T + -^(ea(k) - /i^)(5A^<,(k) 



k,(T 



L ^ /^.,(k,k')5A^.(k)<5A^,,(k'). (1) 
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k,cr:k',o-' 



Here 5A^o-(k) is a small change in the ground-state mo- 
mentum distribution of species a, Ug- is the particle den- 




FIG. 1: A Galilean transformation (dotted curve) and Lorentz 
boost (dashed curve) of a Fermi sphere (solid curve). The 
boosts have been scaled in order to compare the shapes of 
the Fermi sphere after the transformation. The volume of the 
Fermi sphere, i.e., the density of the gas, becomes larger due 
to the Lorentz contraction. 



sity, ecr(k) is the energy of a quasiparticle that for our 
purposes only depends on the magnitude of the vector k, 
and ^a is the chemical potential. In the case of electrons 
the index a specifies the spin of the electron, while for the 
case of quarks it specifics the spin, color and flavor of the 
quark. Hence, the quasiparticle energy and the effective 
interaction between quasiparticles are defined as 
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An important quantity in Fermi-liquid theory is the 
effective mass tti* of a quasiparticle of species a near its 
Fermi surface. It is defined in terms of the group velocity 
V* evaluated at the Fermi momentum k^ 



9e.(k) 



1 -1 



dk 
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such that upon linearizing the quasiparticle energy 
around the Fermi surface we obtain 
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M. 



+ ^(fc-fc<,). 



(3) 



with /i^ = Ecr (fccr) as tlic Fcrml energy of species a. In the 
non-interacting limit, the effective mass reduces to m* = 
Too- and ?Ti* = ea-oikcr) = v'fc2 -I- TO^ for non-relativistic 
and relativistic dispersions, respectively. 

Using the fact that the pressure is an invariant under 
Galilean or Lorentz transformations (see Appendix. IB ip , 
of which the effect can be written as a change in the dis- 
tribution as shown in Fig. [TJ Eq. ([T|) should give zero for 
this particular choice of SNc^ (k) and hence relate e^r (k) to 
/crcr'(k, k'). In other words, by performing an infinites- 
imal Galilean or Lorentz transformation on the Fermi 
sphere, the effective mass can be related to the effective 
interaction at the Fermi surface [l|, [2g . 

To this end, consider an infinitesimal change in the 



Fermi surface from k^ to k„ + 5fco-(k) 



5N^{k) 



K < k < k„ + Sk„{k) Sk„ > 0, 
k„ + Sk„{k) <k<k„ Sk„ < 0, 
otherwise, 



where the change in the Fermi surface solely depends on 
the direction in momentum space denoted by k. The cor- 
responding change in the thermodynamic potential den- 
sity can be written as 



^ ft'fT ~\0 flffj 



dk dk' 
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I dkk^ I dfc'fc'2/<x.'(k,k')-^n,JM. 



T.^J t'^'^^^^T.^^. P^^^f^-^^^^^-^^^^^^^ 



(4) 



r 



In the second line the integrals were expanded to second 
order in 5ka (k) , and fcrcr' {ka- , k^i , 9) was defined as the 
effective interaction at the Fermi surfaces between two 
species, with 9 the angle between the directions of k and 
k'. Next we expand Sk^r in terms of spherical harmonics 
and /o-cr' in terms of Legendre polynomials as 

^fc,(k) = ^<5A:^y,„(k), 

I. in 

fa<7'{ka,k„',0) = ^//,^^'(fc^,fc^/)P;(c0s6'). 
I 

For an imbalanced system, the difference in density ria- 
results in a difference in fia- , which can cause m* to be dif- 
ferent from each other even if their non-interacting bare 
masses nia are the same. For each species, the particle 
density obeys n^ = k^/Qn^. Inserting these expansions 
and using the orthogonality of the spherical harmonics, 
Eq. ^ becomes 
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(5) 



In the balanced case, all Fermi momenta are the same, 
ka = kF, 6k^J^ = Sk^f", then Eq. © reduces to 
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(6) 



where the intensive quantities m* , /i, kp and Fi are the 
same for all species, and n is the total density, i.e., the 
density summed over all species. The (dimensionless) 
Landau Fermi-liquid parameters Fi are defined as 
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kpfl^aa'ikp) 
2^ 



ae(k) 



dk 
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(7) 



Note that Fi is species independent, but the inter-species 
and intra-species interactions can be different, so the sub- 
script a is shown explicitly in fi_aa' but not in Fi. From 
the above it is possible to derive the equation for the ef- 
fective mass, as we first show for the simpler balanced 
case. 

Consider a change in the distribution due to a Lorentz 
transformation as shown in Fig. [T] The shape of the 
Fermi sphere is determined from the equation /i = 
— fc^jU^, where k^ = (e(k),k) is the four-momentum, 
u^ = (7, 7v) is the four- velocity of the Fermi sphere with 
the Lorentz factor 7 = l/Vl — v^ ■ The minus sign in the 
previous expression and in u^u^ = —1 are a consequence 
of our choice of the metric r\^^^ = diag(— 1, 1,1,1). In the 
rest frame, where u^ ~ (liO)j the above reduces to the 
condition £{kp) = ^. Similarly, the shape of the Fermi 
surface can be determined from the Lorentz transforma- 
tion of the momentum to a frame moving with velocity 
-V, 

k ^ k - v(v • k)(l - 7) + e(k)v7, 

which for an infinitesimal Lorentz transformation reduces 
to k — > k + e(k)v and in the non-relativistic limit to 
k — > k -I- mv. Note that the deformation of the Fermi 
sphere due to the Lorentz contraction shows up only in 
higher-order terms of v. 

First, consider the non-relativistic case in which an 
infinitesimal Galilean transformation shifts all momenta 
from k to k -I- mv. This results in a change of the Fermi 
surface according to Skp = rnk ■ v = mvYio{cos9)/\/3, 
such that 6k^p^ = 5iiS,no?nv/^/3. Furthermore, the total 
energy of the system transforms as i? — ?> E + ^Mv^, 
where M is the total mass of the system, which induces 
a change in the chemical potential fi ~ dE/dN -^11 + 
^mv"^. Inserting this into Eq. ([5]) gives the well-known 
Landau effective-mass relation in the balanced case [ll-Q 



TO = 771 



1 



7^1 



(8) 



Secondly, for a relativistic system the momenta are 
shifted by k — )• k-|-e(k)v such that Sk'f^ = SiiSmofiv/\^. 



Since the total energy of the system transforms as E ^i- 
'yE ~ i<^ + \Ev'^ , the change in the chemical potential 
is (5/i — ^iJLv"^ . Inserting these into Eq. (ISl) gives the 
relativistic Landau effective-mass relation 12611 



m* = ^ ( 1 + ^Fi 



(9) 



Note that in the right-hand side Fi also depends on /x, 
because the Fermi momentum kp is a, in general compli- 
cated, function of ji. 

Now for the more general imbalanced case. Since 
the Lorentz transformation is applied uniformly to each 
species, according to the symmetry of the formula, we can 
expect that the contribution from each species in Eq. ([5]) 
should vanish, namely, 



n^fJ-lv^ Sr^ fl,crcr'iK,k^')klkl,fla■^J.a'V ria-flaV 



0. 
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^ a' 

Then we get the general expression for the effective mass. 
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E 



fi.aa'{K,K')N„'{Q) Ha'k„,m* 



^J■a 



(10) 

where A''cr(0) = m*fccr/27r^ is the density of states at the 
Fermi energy of species a. The arguments in /i empha- 
size that, in the imbalanced case, the interaction may 
depend not only on the angle between the momenta but 
also on the absolute value of the Fermi momentum of 
each species. In fact, the same expression of m* can be 
obtained more strictly by considering the addition of a 
single particle to the system and comparing the energy 
increase in two different frames. 

We can introduce the average effective mass, which is 
of practical importance, as m* = '^^m'^/N with N the 
total number of species. Correspondingly, the average 
chemical potential is // = X^cr/^o-/-^- ^^ order to get a 
similar expression as in Eq. 0, we introduce the relative 
weight of each species as x^- = fi^/Nfi such that ^^ Xa- = 
1 . and then generalize the Landau effective- mass relation 
with average values as 



1 



l + xE/i---'(^-'^-')^-'W 



^a' '^a'T^T'o 



(11) 

where the double-sum term, which can be defined again 
as Fi, plays the role of the Landau parameter in the 
balanced case as in Eq. (|9]), therefore we obtain now 



i^l=E/i.--'(^-'^-')^-'(0) 



kam*, 



which, obviously, reduces to Eq. ([7]) as the system be- 
comes balanced. 

In the following sections the effective mass will be de- 
rived from a microscopic RPA calculation. In particular, 
in Sec. IVII it will be determined from the specific heat. 



The specific heat (per volume) for low temperatures can 
be expressed in terms of the effective masses by 



1 fdE 



katn 



V,N 
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dE 9iV^(k) 
aAr^(k) dT 
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V.N 

(12) 



where it was used that d{E/V)/dN„{\i) = dn/dN^(k) + 
ficr = Ecr (k) and for the temperature derivative the Som- 
merfeld expansion was used, c.f. Eq. (|E8[) . This result is 
a trivial generalization of the single species or balanced 
case, since the effective mass already contains the effect 
of interactions, such that the contribution to the heat 
capacity of a single quasiparticle is additive. 

III. RANDOM-PHASE APPROXIMATION 

Our starting point is the partition function 



Z 



VipVi^VA^jyfjVrj exp 






containing the Euclidean QCD Lagrangian density fixed 
in a linear gauge f^^A"^ = 0, 

-Cqcd = E V^/(7m^m + "V - 7oM/)V^/ + \gI,G% 
f 

where ipf is the quark field of fiavor / (the color c and spin 
s degrees of freedom are not shown explicitly) with mass 
•nif and chemical potential /xy, A"^ are the giuon fields, 77° 
are the ghost fields, g is the QCD coupling constant, S, is 
a gauge fixing parameter, f"'^'^ are the fine-structure con- 
stants of the color S'C/(3) group, the covariant derivative 
F>^i = d^ — igt'^A'l^ and the antisymmetric gluon field ten- 
sor G°^ = d^^Al-dyAl+gf^^A^^Al. For large densities, 
the fermionic degrees of freedom become increasingly im- 
portant, such that the large- A^j;- limit could give insight 
into the behavior of quarks at high densities. Therefore, 
consider the large- A^j;^ limit with a fixed rescaled coupling 



-QCD 



-Conn = E ^/ ('^ + "^/ ~ ^^^S - SN^'^ij^f^A"^ ) ijf 
f 

+ \k [d^d, - d'b^, - i/^/.) Al 

1 
2^ 



^gNp-{d^Al~d,A'l)r''^A'>^Al 



\gNp'r'^r^^AlAlAlAl, 



Integrating out the fcrmions and the ghosts gives 



■'-QCD 



Trc./,s In 



G^I[l + gNp^-Gofil^f 



A' 



- Tr, In [d^f^d'^' + gNp-^r^'^A^f^ 
+ \gN~^- {d,Al - d^Al) r'^A^Al 



.^^Ifabcfade^b^^c^d^e^^ 



where the subscripts {c, /, s} exphcitly indicate that the 
trace should also be taken over color, flavor and spin 
space. When we expand the first line in Np, the first- 
order contribution will vanish due to conservation of 
color, i.e., Trc[t°] — 0. Only the second-order expan- 
sion in Np of the first line will give a contribution. Thus 
expansion in powers of Np gives 



''-QCD 



Ti-cj,s ln(-G-;) - Tr, ln(5^/,,<5'^^) 



where we defined the polarization tensor H'^'' from the 
contribution of the quark with fiavor / as 



n°* = ^rooooooo 




'00000000^ 



FIG. 2: The polarization tensor n"^. to the first order. 
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FIG. 3: The thermodynamic potential in the large-A'^F or 
RPA approximation where the double gluon line signifies the 
dressed gluon propagator. 



Sec. lIVI an example is considered to check the expression 
for the Landau effective mass analytically and in Sec. IVII 
the above thermodynamic potential density is calculated 
numerically for the case of a two-flavor imbalanced quark 
system. 



n?:^. = -g'Tr,.AGon,rGon.n 

which is shown diagrammatically in Fig. [5] The polar- 
ization tensor H'i'' is diagonal in color space, as can be 
seen explicitly form Eq. (jEip . Upon integrating out the 
giuons the thermodynamic potential density is 



n{T,{i,f})=-V-'p-HnZ 
1 

, Ng 



NcTrf^sH-Go^) + 



-Tr 



ln(Z?o-^)-21n(a^^)-ln 



N 



^Trln l + i?oEn/ 



(13) 



where Nc is the number of colors, Nq = Nq 



1 



is the number of giuons, and we use the QCD values 
Nc = 3 and Nq = 8 in the following. The first two 
terms arc the ideal Fermi and Bose gas contributions 
to the thermodynamic potential density, where the term 
Trln(9^/^) cancels the two unphysical degrees of freedom 
from Trln(Z3g"^). The last term is the RPA correction 
(ring sum) to the thermodynamic potential density, as 
is shown in Fig [31 Since we are mainly interested in the 
temperature dependence of O and the fact that the above 
definition contains (divergent) zero-temperature contri- 
butions, the T = expression will be subtracted, i.e., 
we consider An{T, {fij}) = n{T, {/i/}) - n{0, {fif}). In 



IV. FERMI-LIQUID PARAMETERS IN WEAK 
COUPLING 

In the limit of weak coupling the interaction contri- 
bution to the thermodynamic potential density is given 
solely by the exchange diagram 



n,, 



Ng 

2V(3 



ETrPo%)- 



Using the definition of the free gluon propagator in the 
Lorentz gauge, this can be written as 



^int -^TJT^Y^ 
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J Sl,S2,S3 — ±1 



(p>q) 



x[7Vf(q)7V^^(p-q)(l-7V;nP)) 

- iVf (p)(l - TVf (q))(l + iVf (P - q))], 

where the Fermi and Bose distributions are defined as 
1 --„ , > 1 



iVf(p) 



p^[s(!o/(p)-pj] ^ Y 



, N%{V) 



3,3[se3(p)] _ 2' 



respectively, the explicit expression for H^.^^ is shown in 
Eq. (|E3|) . the Matsubara sum over the bosonic frequen- 



cies ujq was performed, eo/(p) = Jv^ + mj denotes the 



free quark dispersion while e^(p) = p is the free gluon 
dispersion, and the function J^ was defined as 



j'i^.s^.sAP'^) 



S3 



eo/(p)eo/(q)e^(p-q) 

-eo/(p)eo/(q) + siS2(p-q + 2TO2) 



Sieo/(p) - S2eo/(q) - S3e9(p 



q)' 

(14) 



which will soon be shown to play the role of the inter- 
action between two quarks (si = S2 = 1), two anti- 
quarks (si ~ S2 = —1) or a quark and an antiquark 
(si = — S2 = 1 or si = — S2 — —1) mediated by the emis- 
sion or absorption of a gluon (depending on the sign of 
S3 relative to si and 52)- 



The quasiparticle energy and effective interaction can 
be obtained by varying the thermodynamic potential 
with respect to the Fermi distribution. In the flavor- 
imbalanced case, the distributions depend only on the 
flavor but not on the color and spin indices. To distin- 
guish between the different contributions for each particle 
species it can be checked that the multiplication factor 
No/4: could be written as ^^ ^, ^ ^, t'^c''^'c'c/^^ where s and 
s' are just dummy indices denoting the spin degrees of 
freedom, which are helpful to arrive at the following ex- 
pressions. Note that ilint is independent of the spin in- 
dex, and in fact the spin degrees of freedom has already 
been summed over in the expression for J^. Explicitly we 
find for these quantities 
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ec,/,s(k)-M/ 



sn 



^KfA^) 



/{c,/,4;{c',/%s'}(k,k') 



S' C' S2,S3 = ±1'' ^ ' 

x[(l+iV^^(k-q)(l-7Vf(q))+7V^^(k-q))iVf(q)], (15) 
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E -^i.+,.3(k,k'), 



(16) 
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where Ec' Kc'K'c = ^ and t^,,i^,, = 1/2 - 5,,,l2Nc. 
according to the definitions in Appendix I A II 

In the quasiparticle energy Eq. p^ the first term is the 
noninteracting part, while the second term corresponds 
to the quark self-energy, as is shown in Appendix [G] The 



angular averaged interaction parameters are defined by 
the coefhcients of the Legendre polynomial expansions of 
the effective interaction Eq. (|T6| evaluated on the Fermi 
sphere, by denoting cos 6* = k • k'. 



fl;{cj,s} - E h{cj,s};{c'. 



21 + 1 
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which gives for the zeroth and first interaction parame- 
ters 



/o;/ 
/i;/ 
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8Nc Mg/ 
g^NaNs 3 A^ + 2m) 

8Nc 14] 2k) 
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A2 + ik) 
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A2 + 4k) 



, (17) 



where Ns = 2 is the number of spin degrees of free- 
dom, and a regulatory gluon mass A was introduced, i.e., 
e^ = V^^ -f A^. We see that in this approximation the 
interaction is color- and spin- independent. For the sake 



of simplicity, therefore, we from now on omit the unnec- 
essary color and spin indices, and use kj and ^0/ for the 
Fermi momentum and the non-interacting Fermi energy 
of the quark with flavor /. The effective mass of the bal- 
anced case, Eq. ^, can be obtained from the interaction 
parameters Eq. (jT?)) . and Eq. ^, by expanding to the 
first order in the coupling constant 
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fi kpm* 



'^ 3 27r2 
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M^ 



67r2 



/i-f-0(/), (18) 



where it was used that fi = fio + 0{g^) as follows from 
Eq. ([T5|). For the imbalanced case, a similar result can 
be obtained for each species according to Eq. ((TO)) . 



^/ = Mi- 



^^ffl■.Akf)Nf{o) 



A^/+^/i;/(fc/)' (19) 



where the Kronccker delta Sff in the interaction simph- 
fies the expression considerably, such that we obtain the 
similar expression as in Eq. ()18|) . 

As a special case, we can use the above results to dis- 
cuss a two-flavor imbalanced system, where the average 
effective mass and chemical potential can be defined as 
m* = {m*^ -f m*_)/2 and ^ ~ {fi-^- + fJ--)/'2 with the 
subscripts -I- and — for the majority and minority fla- 
vors, respectively. The imbalance can be quantified as 
h = (/x_|_ — /i_)/(/i+ + A*-). For small imbalance, the 
factors k±fi.± can be expanded in h with respect to 
the value in the balanced case, namely k±fi.±(k±) « 
kFfi{kF){i ± Ai±/i -|- A2±h^). Note that, the linear co- 
efficients of the expansion can, in general, be different for 
each flavor. We obtain the expression of m* for small h 



m 



kpi^^fiikp 



(2 



127r2 

2Ai+ 



Ai_)/i 



-[2 + (An 

2Ai_ + A2+ + A2-)/i2 + 0(/i')], 



where the h^ term is kept since the linear term may van- 
ish for the symmetric case, namely Ai+ = Ai_, and 
we have replaced ^o with // because the difference is in 
higher order of g. Compared with Eq. pTj) . the effective 
Landau parameter for the average effective mass is 



+ (2 + 2Ai+ 



(Al4 

2Ai_ 



Ai_)/i 

- A2+ + A2- 



0{h% 



Similarly, we can obtain the effective mass difference 
Am* = ("T-+ — 7Til)/2 as a function of h 



A,„. =^ih 



kpfJ-'^fiikp) 



[(4- 



127r2 

(2Ai+-2Ai_+A2+ 



Ai+-fAi_)/i 



It is quite natural to find A^* (x h in the leading order. 
The consistency of the above results will now be ver- 
ified to the lowest order in the coupling constant using 



the quark self-energy. Starting from the dispersion and 
self-energy of Eq. (|15p , the effective mass in Eq. ^TE\\ or 
Eq. (|19p can also be derived in a different manner. The 
quasiparticle pole of the quark propagator is 

e/(p) = eo/(p) + S+/^,(eo/(p) - M,P), 

where E+ is the renormalized positive-energy projected 
self-energy. The effective mass corresponding to this pole 
is most easily defined from the Fermi velocity of the 
quasiparticle 



V _ «/ 



M/ 



"f'f 



1 



^e/(p) 
e/(p) V 9p 



J p=kf 



Inserting the definition of the quasiparticle pole into the 
above gives, using deof{p)/dp = p/eo/(p), 



=M/^/ 



eo/(p) 



9eo/(p) 



=M/ 



dp 

deQf{p)J:+jJeof{p) - pt/,p) 
dp 

(9eo/(p)i;+f,(eo/(p)-M/,P) 



0{g' 



Jp=fc/ 



9eof{p) 



+ 0{g'). 



^p=kf 



(20) 



Using that the explicit form of the self-energy in the zero- 
temperature limit can be written as (c.f. Eq. (|G2p ) 



^c;/,.(£o/(p)-M/,p) 



d^g 

J^f{cJ^s}:{c'J'^s'}{P,q.)N+{q), 



E 

c'J',s' 



and the relation Eq. (jBip between the vector derivatives 
of the effective interaction, the derivative in Eq. (|20p can 
be rewritten as 



J 



dp deaf{p)T.+fJeof{p) - M/, P) 



5eo/(p) dp 



f/l?'^'"'«»*(-'4 



Jp=fc/ 






E 

c',f',s' 



d'9 . g£0/(p)/{c,/,s};{c'./',s'}(p,q) ,,+ ,„, 

3P 7^ ^V' W 



(27r) 



3/ i^P-^ E f{cJ,s};{c'J',s'}{p,Cl) 



c'J'.s' 



p=kf 



dp 

kf^J■lf fl;{cj,s}{kf) 

27r2 3 



p=kf 



where partial integration was used and Sff'dqNt,{q) = microscopic diagrammatic calculation. 
—Sffiq5{q~ kf) at zero temperature. Therefore Eq. ([20)) 
is identical to Eq. (fT9|). which shows that the phenomi- 
nological Landau argument is indeed consistent with the 



V. DRESSED GLUON PROPAGATOR 

Before we start the calculation of the temperature de- 
pendence of the large- A'^p thermodynamic potential, it is 
useful to examine the dressed gluon propagator, which 
is a crucial ingredient of the RPA theory. The most 
physical gauge to study the propagator is the Coulomb 
gauge, since its form results from considering linear re- 
sponse [12, \lM ■ This is also by far the most used gauge in 
condensed-matter theory. In this gauge the gluon prop- 
agator is |12 | 



D^.iQ) 






Q Sfj^oSiyQ 



02 + G{Q) q'- Q 



F{Q) ' q^ Q2 



g" 



Q' 



where Pj^ is the three-dimensional transverse projector 
defined in Appendix IA21 F and G are related to the 
longitudinal and transverse projections of 11^^ including 
the contributions of various species as well as the vacuum. 
More details can be found in Appendix |E1 

Consider the spectral functions of the transverse and 
longitudinal propagator, which are defined by 



«^'^(a.,q) 



-^ 



[D^'\u: + iO,q)] 



In general the spectral function depends on the gauge, 
however, the positions of the poles are gauge indepen- 
dent, and the Coulomb gauge has the additional property 
p^'^(cj, q) > for w > as required of a physical spectral 
function. For the balanced case, the spectral functions 
are shown in Fig. U] for several values of ^ and in Fig. [5] 
for several values of T. The latter includes a small-T 
correction to F and G, see Eq. (|E9p [45| . The transverse 
and longitudinal plasmon modes and the large contribu- 
tion due to the decay of the gluon into the particle-hole 
continuum {0 < oj < q) are clearly visible. Note that 
massless quarks are used in the limit of high density, c.f. 
Appendix [E] However, a small but finite quark mass is 
necessary to renormalize the real part of the polariza- 
tion tensor of the vacuum, as shown in Eq. (jE7p . Since 
the vacuum contribution plays no significant role in the 
following calculation, we simply take this nonzero quark 
mass m in the vacuum term the same for different flavors. 

For the two-flavor imbalanced case, it is clear from 
Eq. p3p that in the leading-order correction each flavor 
contributes separately to the polarization tensor. Be- 
cause of the above setting, there is no mass imbalance in 
the present system, therefore it is enough to consider only 
positive h due to the symmetry. Following the notation 
in Sec. IIVI all the above results can be easily generalized 
to such an imbalanced system by using ii{l + h) to re- 
place /x+, and fj,{l — h) for fi-. As /i — ?> 0, the system 
reduces to the balanced case. Now the low temperature 
condition requires T <IC /x(l ± h), such that h can not 
be too close to 1, namely the extremely imbalanced case. 
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FIG. 4: (color online). The spectral function of the transverse 
(upper panel) and longitudinal (lower panel) gluon propaga- 
tors as a function of frequency oj for several values of /i in the 
balanced case. Shown here for T = Q,m = 1, q = 20, g = 1/2, 
Nf = 2 and for different curves /i = 1,20,50,100,200,300 
from red to purple. 



The spectral function for a two-flavor imbalanced system 
is shown for various h in Fig. |6l Comparing with Fig. [5l 
we see that increasing h has a similar effect as increasing 
T. 

The dispersion relations of the modes can be found by 
solving 



^ 



Dt\{^,(i) 



= 0. 



(21) 



In the transverse case there is a single solution ^'^ (q) > q, 
however, in the longitudinal case there are two solutions 
^f (g) < q < £.2{q), as shown in Fig. [3 Note that ^f (g) 
is not a real propagating mode because the imaginary 
part in the region < to < q is large due to particle-hole 
creation processes, as explained in Appendix [E] The two 
plasmon modes £,^{q) and £,2{q) approach the so-called 
plasma frequency Wpi as g — >■ 0. The plasma frequency 
of the transverse and longitudinal mode can be found by 
expanding the inverse propagators for small q and small 
w, which in the zero-temperature limit gives for both 
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FIG. 5; (color online). The spectral functions of the trans- 
verse (upper panel) and longitudinal (lower panel) giuons 
as functions of u) for various T in the balanced case, with 
/x = 200, m = 1, q = 20, 5 = 1/2, Np = 2 and for differ- 
ent curves T = 1, 40, 80, 120, 160, 200 from red to purple. (In 
principle, T should be much smaller than fi, but here we show 
also large values of T to demonstrate the curves more clearly.) 



cases 

1 T 1 ^ n — ^ =0, 

3 ^2 V 4/x2 4^2 j 

whose solution for small coupling constant yields Wpi ~ 
\/2/2>mg [12| . Here rUg is the gluon thermal mass whose 
expression is obtained in the hard dense and hard thermal 
loop approximation as shown in Eq. (|E10[) . In the limit 
q ^ rag ^ gfi all solutions reduce tooj = q, sec Eq. (JEIII) . 

Furthermore, for the limit u),q ^ fi the vacuum be- 
comes increasingly dominant such that Eq. ([?!]) has 
a zero at large Q called the Landau pole, Q^ = 
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FIG. 6: (color online). The spectral functions of the trans- 
verse (upper panel) and longitudinal (lower panel) giuons as 
functions of cu for various h, with T = 0, m — 0.001, q = 0.1, 
/i = 1, g = 1/2 and for different curves /i = 0, 0.2, 0.4, 0.6, 0.8 
from red to blue. 




expff -I 
however 



247r'' 



m 



2 — 



^ M 



plays no role for our purpose as we are inter- 
ested in the low-temperature behavior of the theory that 
is hardly influenced by the high-energy behavior of the 
gluon propagator. The dispersion relations for the case 



c.f. Eq. dEllD, which, pjQ 7. ^pj^g solutions of Eq. ([2T]) at T = /i = 0, with the dot- 
dashed curve for (^^ , the dashed curve for ^f , and dotted curve 
for ^i . The solid line is a guide to the eye and corresponds 
to the dispersion of undressed giuons, oj = q. 
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with nonzero T ot his similar but just with a httlc higher 
TOpj. In fact, as a generahzation of Eq. (|E10p . the thermal 
mass of a two-flavor imbalanced quark system reads 



9 

47r2 



^l^{l + hf + ^i^{l-hf + 



27r2T2 



where h thus plays the same role as ttT/vS/^. 



VI. THE FULL LARGE-TVf THERMODYNAMIC 
POTENTIAL 



The full ideal-gas contribution to the thermodynamie 
potential density with the zero-temperature contribution 
subtracted is (c.f. Appendix [F| 



Afin = -Nr 



Nf 



180 






'^'f 



Nc 



rp4 

'45^ 



The first part is the contribution of an ideal massless 
Fermi gas, while the second part is the Stefan-Boltzmann 
law of an ideal Bose gas. The RPA correction to il, the 
last term in Eq. (|13p . can be written as 



^RPA(r, {/i/}) 



No 

2VP 



^Trln(l + i:'on) 



^El-Det(l + ^^ + ^ 



l^n ,(J 



l|i:;»0.^)0.|). 



where the Lorentz gauge is used, while for the Coulomb 
gauge the second identity is not valid but the last result 
is still the same, which is a consequence of the gauge in- 
variance of the thermodynamic potential density. Using 
contour deformations to carry out the Matsubara sum, 
as shown in Appendix [Dl we obtain 



n 



RPA 



miM) 



No 
2tt 



dqduj 



'A7Tq^[2NBiio) + l] 
{3[ln^(c.+ ,q)]+23[lnG(c.+ ,q)]}, 



(2^)3 



where F(w+,q) = 1-1- F(a;+,q)/(— w^ -|- q"^), and sim- 
ilar for G by replacing F with G. The temperature 
dependence can be obtained by subtracting the zero- 
temperature contribution. Since F and G contain correc- 
tions of order T^, the leading-order correction in AfjRPA 



comes from two parts, 
Ar2R,pA(r, {^/}) = ^ / dqdujq'^ 



Nb{uj) 



c 



arctan |||^ + 7re[-5R(/'°)]sgn[5(F)0] 
-2arctan^^ + 27re[-3fJ(G'°)]sgn[3(G°)] 



5R(G0) 



- [3(ln F' ) - 3(ln F") + 23(ln G' ) - 23(ln G")] 



where the superscript T or means the corresponding 
terms are taken at nonzero T or T = 0. We will refer 
to the first part as the Nb term and the second as the 
non-NB term. The leading correction from the non-iVs 
term can be shown to be proportional to T^ and is not 
of great interest in our study, since we will concentrate 
on the anomalous and dominant T dependence, which 
is a consequence of the Nb term. In the Nb term, the 
arctangent terms can be interpreted as contributions due 
to production and decay of thermal gluons, because of 
their dependence on the imaginary part of the gluon self- 
energies F and G, while the theta function terms are 
interpreted as a correction to the ideal gas law due to 
thermal plasmon modes. 

The frequency integral over the theta function can be 
performed explicitly. In both the transverse and longitu- 
dinal case the sign of the imaginary part is positive when 
the real part is negative, such that after the frequency 
integration the result of the integral is proportional to 



rp4 



X In 



1 _ e-/3«2(/3-':r) 

1 _ p-/3«i(/3-i2:) 



dx, 



where x = q/T and S,i < i^ < S,2 signifies the region where 
'^F,G < 0. In the limit of low temperature this integral 
will go to a constant. The theta- function contribution 
can thus be neglected since it is of higher order in the 
temperature than is of interest to us here. 

Next we perform the integrals over the arctangents, 
whose structure at T = /i = is shown in Fig. [Sj Note 
that, for the study of the leading-order T corrections, it 
is not necessary to include the T^ term in F and G since 
the integral with Nb(j-^) at low T is already in the order 
of r^. The dominant contribution for small temperatures 
(T <^ ii) comes from the frequency integration over the 
domain w e [0, q] for the case q < 2/i, which is due to 
particle-hole creation. For the two-flavor balanced case 
{Np = 2, h = 0), it was found numerically that the 
integral in the limit of small temperatures is 



^2 



(cf-4lng2) 



for the electric (longitudinal) gluons and 
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FIG. 8: (color online). The structure of a,Tcta,n{Q f/^Hf) (up- 
per panel) and arctan(Q'G/3?G) (lower panel) with T = 0, 
m — O.OOlfi, h — 0, g — 1/2, and for different curves, 
q//i = 0.8 (red), 1.6 (green), and 2.4 (blue), respectively. 
The contribution for max(0, q — 2/i) < uj < q is solely due to 
particle-hole contributions, while the contribution g < a; < oo 
is due to particle-antiparticle processes (a combination of fi- 
nite density and vacuum processes). The discontinuities at 
(jj = q correspond to the plasmon modes shown in Fig. [7] 
The constant tails are due to the vacuum contribution, which, 
however, will not cause divergence because of the Bose distri- 
bution function Nb{l^)- 



for the magnetic (transverse) gluons, with 



0.16, 



0.42, and 



~ 0.48, 
0.056. These results are 



quite close to the analytic results for the leading-order 
correction terms at small g and T obtained in Ref. [35| . 
where cf = [ln(47r2) - l]/6 w 0.4460, ^ = 1/6 w 0.1667, 

cf = [jE - 6C'(2)/7r2 + 3/2 + ln(327r)]/18 w 0.4032, and 
c'^ = 1/18 « 0.05556. 



For the imbalanced case, the above result for the trans- 
verse gluons can be generalized to 



s=±l 



In- 



9 



nil + sh) 



(22) 

where the contributions from both flavors with chemi- 
cal potentials fi{l ± h) are additive. The longitudinal 
part needs some further discussion. Since the static long- 
wavclcngth longitudinal modes are screened (c.f. Ap- 
pendix [E| , the two chemical potentials contributing to 
^F can not be separated even in the lowest-order term. 
This is different from the result obtained in weak cou- 
pling, as shown in Eq. (|16p , where no interaction between 
the two flavors is involved. Therefore it is not surprising 
to find that a simple generalization of the balanced case, 
as Y.s=±i g^fJ''^{l + sh)^T^{c^~c^ In.g2)/27r2, does not fit 
well with the numerical results. To obtain a reasonable 
ansatz for the longitudinal part, wc integrate the corre- 
sponding arctangent term of the imbalanced case in the 
low temperature limit up to 0{g^) to obtain 



■E 

s=±l 



5V(l + sM'r2 



1 47r2(l 
6^" 7^ 



shf 



h^) 



where wc sec that the factor (l + /i^) = | S.s=±i(l' 



(23) 
'Shy 



is a mixture effect of the two flavors. This expression fits 
the numerical results very well with an error of only 3.5%. 
Furthermore, just based on the numerical data, we find 
another ansatz which fits the results even better 



Ij: 



5V(1+ 5/^)2^2 



=±1 



lu- 



ff 



1 + sh 



(24) 



where the mixture effect is shown implicitly in the log- 
arithm, since the denominator becomes dimensionless 
by canceling with the average chemical potential /i = 
^g^^j^ /i(l + sh)/2 in the numerator. 

The specific hea t (p er volume) at fixed volume and 
particle number is |29| 



Cv 



dT 



= T 



dS_\ 
df) 



M/ 



/ 



{dnf/d^if)T 



where the entropy density S — {dfl/dT)v and parti- 
cle number density nf — —(dn/dfif)v- In the low- 
temperature limit the second term can be neglected, 
therefore the specific heat can be obtained, using 
Eqs. (1251) and (EII), as 



Cv-C^^ -T 






M/ 



■E 

s=±l 



2gV'(l + s/i)'T 



1, 47r2(l 
- in — Try — 
6 5^(1 



shy 



h^) 



2cf + 3cl 



2c2 In 



g'T 



Ai(l + s/i) 



(25) 
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'^\^ of an 



where the specific heat C^ = -Tid^Afla/dT^) 
ideal gas has been subtracted. 

The effective mass and the first Landau parameter can 
be determined in the high-density limit by comparing 
Eq. (O and Eq. ^. 



^J■± 



g^ti{l±h) 



In 



9^1 + h^) 



^12cf + 18c^ + 12c^ln 



1 



g'T 



/i(l±/i) 



(26) 



where we used k± — fi± = fi{l ± h) in the massless 
limit at high density. The first term in the effective 
mass is due to the ideal-gas specific heat. However, as 
pointed out before, the factor (1 + /i^) shows the mixing 
of the two flavors, unlike the weak-coupling results ob- 
tained in Sec. IIVI e.g., Eq. p9)) . This is because, with 
the RPA correction, we incorporate the sum of an infi- 
nite chain of gluon self-energies, which incorporates in- 
teractions between different quark flavors. Since the lon- 
gitudinal gluon is screened in the static long-wavelength 
limit, the mixing is even present in the low-T and small-5 
limit. The numerical fit in Eq 1^^ also provides another 
expression for the effective mass. 



TO_i 



^J■± 



.g2^(l ± h) 



6cf - 6c^ In 



9 



-12c?' + 18c^ + 12c^ln 



l±h 
g'T 1 


fi{i±h) 



(27) 



As expected, all the above results return to the balanced 
case as /i — > 0. As mentioned previously, the contribution 
from the non-Ns term is not included, which acts as a 
constant shift on the cf and cf factors. In Fig. [9] the 
change of effective mass due to the interaction is given as 
function of imbalance at various temperatures. Finally, 
we emphasize again that h should not be too close to 1 
even though the divergence from ln(l — h) is suppressed 
by the prefactor (1 — h), because in such an extremely 
imbalanced case the condition T <C /i- is not satisfied 
for the minority flavor. 



VII. SUMMARY AND DISCUSSION 

We have calculated the thermodynamic potential per- 
turbatively in the large- A''^ limit and the effective mass 
of the quarks is determined by using Fermi-liquid the- 
ory for an imbalanced cold dense quark system. The 
temperature dependence is obtained by using the gluon 
self-energy, from which the contributions from transverse 
and longitudinal gluons are explicitly shown. For the 
two-flavor imbalanced quark system, the effective mass 
is obtained both analytically within the weak-coupling 
limit, and numerically within the RPA approximation. 
We find that, in contrast to the weak-coupling result, 
where the effective mass of each flavor is independent of 
each other due to the lack of an inter-flavor interaction, 




FIG. 9: (color online). The RPA correction to tlie effective 
mass as function of h, with the red curves for the majority 
flavor (-I-) and the blue curves for the minority flavor (— ). The 
solid curves correspond to the expression in Eq. (|26p while 
the dashed curves to Eq. (|27|l . and their difference is not too 
much, g — 1/2 and for each group of curves the temperature is 
r//i — 10^", 10^* and 10^^ from bottom to top, respectively. 



the effective mass obtained from the RPA calculation de- 
pends on the chemical potentials of both flavors, which 
is a consequence of the static screening of the electric 
gluons at long wavelength. 

From the RPA results, the logarithmic dependence 
on temperature of the specific heat and effective mass 
signals a breakdown of Fermi-liquid theory at zero- 
temperature. Non-Fermi-liquid behavior, due to un- 
screened long-range magnetic interactions, was already 
discussed several decades ago for the case of the electron 
gas [30, ISll- In the large- A^f limit the QCD thermody- 
namic potential is essentially the same as that of QED, 
apart from group theory factors. The non-abelian effects 
of QCD only show up if gluon self-interaction corrections 
are included. The logarithmic behavior in the balanced 
case has previously been seen in analytic and numerical 
calculations of the QCD thermodynamic potential and 
specific heat using the large- Aj^ limit, dimensional reduc- 
tion, and hard-dense loop QCD perturbation theory [32- 
[371 . We expect that such a logarithmic dependence shall 
have important effects on the imbalanced QCD phase 
diagram at low temperatures and directly influence the 
properties of quark matter in the core of neutron stars. 
To better understand these effects, it is necessary to go 
further, such as including the gluon self-interaction and 
extending our discussion to the three-flavor imbalanced 
case. Wc hope progress along these directions will be 
achieved in the near future. 
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Appendix A: Conventions 
1. Euclidean conventions 

The four-momentuni vectors in Euclidean space-time 
will be written with capital letters Q = (i(iw„), q), while 
for the Wick- rotated case (iw„ — >■ w -I- iO) the roman capi- 
tal letters Q+ = (i(w-l-iO), q) are used. Three momentum 
vectors are written in bold face q and its length as q. And 
w ± lO arc sometimes written as (jJ± for short. 

The Euclidean gamma matrices in the standard repre- 
sentation are 



aj.6l 1 Siah 



70 



1 
-1 



li = -I 



a, 
-a. 



75 = 7o7i7273 



1 

1 



All the above entries are 2x2 matrices and di are the 
Pauli spin matrices 



1 
0-1 = I ^ Q ) , cr2 



-i 
i 



c^s 



1 
-1 



The gamma matrices obey the following relations 
{-f,„ji,} = 2Sf,^, {75,7^} = 0, 7I = 1. 
The charge conjugation matrix is 

(72 



C = 7o72 = -i 



0-2 



which satisfies 



C7JC' = 7^, C^ = C'-' = -C. 

The eight generators of the fundamental representation 
of SU(3) are taken to be 



1 



1 

1 




1 



t^ = 



1 



i- = - -1 
^ * 



1 



-i 
i 




1 



r = - 
2*100 





t^ = ^ I 

i 



1 







i° = - 1 
2*010 




2^3 




and have been normalized according to Tr[t°i ] 
In general, the product of the generators is 

Nc - 1 Nc + 1 
^ij ^ki = ~T7rT — ('^*J ^^i + ^«' ^^] ) 7Tt — ('^y' ^^i ^ SaSkj ) . 



4iVc 



^Nc 



Furthermore, one will frequently encounter the following 
group-theory factors 



abc rabd 



r^f 



Nc5' 



N'c-l^ No, 



^a^i3- 2Nc ''" 2Nc '" 

in the above expressions Nc is the number of colors and 
Ng the number of gluons. 

Fourier transforms are normalized as 



^iPf^Xfj, 



where a;„ = 7r(2n + l)//3 are the fermionic Matsubara fre- 
quencies, V = /3{2tt)^ is the imaginary time phase space 
volume and p^ ~ (i(iw„),p) and x^ = (t, x). 
The Dirac equation in Euclidean space is 

{0 + m)iix) =0, 

where $ = 7^9^- In momentum space the Dirac equation 
reads 

{ip + m)ip{p) — 0. 

The eigenvalues of the matrix ip are ±m, since {if)'^ = 
m?. The eigenspinors corresponding to these 



-P 



eigenvalues are 

ipUs{p) = -?71Us(p), 

ipVsip) = mvsip). 

Note that ip is not a hermitian matrix and that Ws(— p) 
satisfies the same equation as Us(p)- However, they can 
also be viewed as the the eigenspinors of a hermitian 
matrix 

Joiip ■ J + m)us{p) = epUs(p), 
7o(ip ■ J + m)vs{-p) = -epVs{-p). 

The positive and negative energy eigenspinors have the 
following form 



Usip) = 



W.s(-P) = 



2e 



6 



P \ Hp+m 
— ^^^^ 



2e 



p 



P I tp+m 



(Al) 
(A2) 
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where ^-f = (IjO)"^ and ^j. = (0, 1)"^. They are orthonor- 
mal in the sense that 

w1(pX(p) = ^'l(-p)«s(-p) = Sss', 
^t(-p)<(p)=«t(_p),/(p)^0, 

and satisfy the completeness relation 

Y,[Us{pWsiP) + Vsi-pHi-p)] = 1. 
s 

Using the above eigenspinors the positive and negative 
energy projectors can be defined as 



s 

p^(p)^^t,,(-pKt(_p) = 



f/ X _ £p7o -ij-p + m 



2e„ 



70, 



ep7o + ^7 • P - "^ 
2e„ 



70, 



which in the massless case reduce to 

1 + sijol ■ P 



P|(p,m = 0) 



2 



The helicity projection operators project the spin along 
the momentum of the particle and read 



^^(P) = 



1 + ,sS • p 1 + si757o7 • P 



a 



where S = 175707 == ( n ^ I is the spin operator. Note 
that[P|(p)7o,P^'(p)]=0. 

2. The Green's functions 

The Green's functions of quarks and gluons in Eu- 
clidean space are defined by 

G(r,x;r',x') = -(V(T,x)^(r',x')), 
I?^,(r,x;T',x') - (A^(t,x)A,(t',x')). 

In momentum space, the free quark propagator is 

{iujn + m)7o — ij ■ P + m if — m 



Goip) 



{iiOn + IJ-)^ — p^ — m? P^ + 771^ 



In the Lorentz gauge [d^Af^ = 0), the momentum-space 
free gluon propagator is 



Do,tii'{q) 



1 



^,..-(1-0- 



Q^ 



n2 '^7)2 n2 ' 



while in the Coulomb gauge {diAi —Q) it is [T2[ 



-Di 



0,/ii' 



^^1/ Q <5^o^f ^ Q Q^Qy 



g2 g2 



e- 



where ^ is a gauge-fixing parameter and the projectors 
can be written as 



P — (5 - 



pi-s.^,- 






pT _ pT _ PI 



pL _ p _ pT 



Note that the free gluon propagator in the Lorentz gauge 
satisfies Q^Dq^^i, = ^Qi^/Q^, such that in the Landau 
gauge (^ = 0) it is purely four-momentum transverse, 
i.e., QpP'o,^^ = 0. The free gluon propagator in the 
Coulomb gauge satisfies q,;Z3o,w = ^q^Qi^/Q^, such that 
in the Landau gauge it is three-momentum transverse. 
These propagators arc diagonal in color space. 



Appendix B: Lorentz transformation properties 

In this section the Lorentz transformation properties of 
some quantities are summarized, such as the thermody- 
namic potential, the distribution function and the effec- 
tive interaction. The Lorentz transformation to a frame 
moving with velocity v is 



/-7 



A'^''(v) 



-7V 



7v 



<5., + ^(7-l) 



V 



\ 



/ 



such that the four- momentum P'' = (e(p), p) of a parti- 
cle transforms as 



<P) 
P 



7(e(p)- v-p) 
p + v(p-v)(7-l)-7e(p)v 



Not to be confused with the Dirac matrices 7^, the 7 
used in this section is the Lorentz factor 7 = l/Vl — v'^. 



1. Lorentz invariance of the thermodynamic 
potential density 

The invariance of the thermodynamic potential density 
under Lorentz transformations can be shown using the 
stress-energy tensor. Consider the change in the thermo- 
dynamic potential density under a Lorentz transforma- 
tion from the rest frame to a frame moving with velocity 
V. Since an interacting gas of quarks in the rest frame is 
specified only by an energy density p and a pressure p, 
the stress-energy tensor is diagonal and of the form 

T^^ =diag(p,p,p,p)^^, 

such that under a Lorentz transformation the thermody- 
namic potential density fl — p ~ pn — —p changes as 

6Q = SToo — 5{pn) = "f'^v^{p +p) — 5[pn) 
= pn^'^v'^ — 6{pn). 
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The transformation of the chemical potential (/i = 
dE/dN) follows from the Lorentz boosted total energy 
of the system E — > ^{E — w-V) = 7_E, where P = is the 
total momentum of the system in the rest frame, which 
results in ^ — >■ 7^. The transformation of the density is 
due to a Lorentz contraction in the volume n — > 771. Thus 
the total change in ^n is 5{iin) = (7^ — l)/in = 7^w^/in, 
such that under a Lorentz transformation 6fl — 0. For an 
alternative derivation which uses that the pressure trans- 
forms the same way as a force per area, see Ref. |38| . 



Note that the distribution function iV(p) is a Lorentz 
invariant, which can be easily derived from the fact 
that the number of particles in a volume d^xd^p of 
phase space is invariant under Lorentz transformations 
j3g |. i.e., N{p) = -/V(p) where the tilde signifies the 
Lorentz transformed quantity. Subsequently, it is pos- 
sible to derive a transformation law for /(p,p') by us- 
ing that the distribution, the thermodynamic potential 
density and dp/eo(p) are Lorentz invariant. It follows 
that eoa{p)£Qa'{p')faa'{p,p') should bc Lorentz invari- 
ant, giving 



2. Transformation properties of /(p,p') 

Consider the lowest-order correction to the free ther- 
modynamic potential density 

f^int = iX;/^^^/-'(P'P')^-(P)^-'(P')• 



eOa(p)eOa'(p')/<T<T'(P,p') = eOa{p)^Oa'{p')faa'iP,p'), 



where eoo-(p) = \/P + rn 

inv, using eo(p) = eo{p)-v-p + 0{v'^), p 

0{v'^), and p/eo(p) = aeo(p)/9p. 



Expand / to the lowest order 
p-eo(p)v-f 



/<T<T'(P,P') =/<T<T'(P,P') 



<^Oa{p)_ 



eo<T'(p') 



+ Oiv^) 



= /a<T'(p,P ) - eo^rlPJV TZ. eOa'(P )V TT-, ./<T<T'(P,P ) 



dp 



dp' 



--Lcr' (P, P') - V ■ 



5eo<T(p)/a<7'(P.P0 , 3eoCT'(p')/<7<7'(p,p') 



dp 



dp' 



+ 0{v' 



deoaip) 9eo<T'(p') 

V • — V 



9p 



dp' 



+ 0{v^) 



r 



If it is assumed that to the lowest order the interac- 
tion does not depend on any distribution functions, i.e., 
faa'{p,p') = fcra'{p,p'), the abovc implies 



diO<T{p)faa'iP,P') 

dp 



deoa,{p')f^a>{p,p') 

dp' 



(Bl) 



Note that the above derivation is similar to that given in 
Ref. M. 



Appendix C: From Minko\vski to Euclidian space 

In this section it is summarized how to turn the 
Minkowski quantum field theory (QFT) of quantum 
chromodynamics into a Euclidean statistical field theory 
(SET) suitable for studying the dynamical properties of 
a many-particle system. The starting point is the gauge- 
fixed path integral for QCD 



^A^^ip^^p^fj^r] expli / Jf( 



'QCD'^ 



(CI) 



where A^ are the gluon fields, ip and ip = iip'^j^ the 
quark fields, rj and fj = irj^^^ the ghost fields and the 
QCD Lagrangian in Minkowski space fixed in a linear 



gauge (/^A^ = 0) is given by 

f 

-lA^^{d,d,-d'v,.)A: 

- fj-d^d^Y - 9r'"'\fi''d^^if)Al 



- \g{d^Al - d^A>:,)r'^i^,,i^,pAlA'', 

- \9^r''T''ii^.v.pA>:,A:A'lA'i + hrAi)\ 



2C 



(C2) 



where ■(/'/ is the quark field with flavor /, and the sum- 
mation over color and spin indices are shown implicitly. 
The metric was chosen to be 77^1/ = diag (—1,1,1,1) and 
the Minkowski gamma matrices in the standard repre- 
sentation are 



7 = — z 



1 
-1 



7^ = ^70-^1-^2^3 



7 = — i 



1 

1 



a, 

-<J^ 



(C3) 



The gamma matrices satisfy 



5\2 



{7^7''} = 2,7^^ {7^7n = 0, (7') 
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By performing a Wick rotation the above quantum 
field theory can be transformed into a statistical field 
theory. A Wick rotation amounts to taking an analytic 
continuation from real time to imaginary time (t = — ir), 



which turns the Minkowski metric ds 



-dt^ + dx into 



the Euclidean metric ds = dr^ + da;^. To do this con- 
sistently the zeroth component of all four-vectors need 
to change accordingly. The procedure is most easily un- 
derstood by considering the length of the position four- 
vector 






^w.^ = (t,x)f -;i)(i 



(it,x)l( '^ ]=x'^S^,x^, 



..-.. = (-t,x)(-i;)(-^^) 



(zt,x)I( 'j )^x^6,,xl 



where it is seen that the minus sign of the Minkowski 
metric is absorbed in the definition of the Euclidean 
four- vectors x^ = («i,x) = (t, x). In Euclidean space 
no distinction is made between upper and lower indices. 
A simple way to obtain the Euclidean form of a vec- 
tor is to multiply the contravariant vector by the matrix 
diag (i, 1, 1, 1) or the covariant vector by diag {—i, 1, 1, 1) 
and set t = —it. Note that the spatial components do not 
change. For example, the Euclidean position four- vector, 
the four-divergence and the zeroth gamma matrix are in 
terms of their Minkowski definitions Eq. (|C3[) 

2^f = («a^°,x) = (-za;o,x) = (if,x) = (r,x). 



af^ = (zaO,V)-(-z9o,V) 



-4'^)=(^-^)' 



7o^ = ^7" 



-«7o 



1 
-1 



Generalizations to tensors is straightforward and fol- 
lows for instance from the example A^j^i, = a^au. Using 
the above procedure to find the Euclidean versions of 
all tensors, the partition function is easily found from 
Eqs. (JC1[) and (|C2[) by setting the tensors to their Eu- 
clidean versions, taking -q^i^ — >■ (5^,^ and t — > —it. The 
partition function is 



Z 



3)A^^2)-^SJiI}3it\3it\ e~ 



where the Euclidean action is defined as 



S"" = / ^QCDdrdx, 



with the Euclidean Lagrangian 



^QCD = II[i'/(^ + "V - 7o/x/)V'/ - igi'n^^i'i^fK] 






+ ^g{d^A:-d,Al)r'^AlAt 



2 rabc rade a b 



g'TT^AlAlAlAl 



2C 



[d.Alf 



In the above the tensors arc all Euclidean but the in- 
dex E has been dropped for convenience, the conjugate 
fields are now defined as ^^^ = iptJo ^-nd f] — r/^^^ and 
the chemical potential has been added as the Lagrange 
multiplier of the density iplipf. Additionally, to com- 
plete the connection between QFT and SET, the time 
integration domain is changed to t G [0,/3], where f3 
is the inverse temperature T of the system. Due to 
the definition of the partition function as a trace over 
all states, the bosonic (fermionic) fields are required to 
obey symmetric (anti-symmetric) boundary conditions, 
namely i^fir = 0, x) = ±V'/(t = /3, x). 



Appendix D: Nonzero temperature calculations 

At nonzero temperature one needs to calculate Mat- 
subara summations which usually can be done using con- 
tour integration. In the following an expression is de- 
rived for such summations and the interpretation of the 
result is discussed. Consider to this end the sum over 
bosonic Matsubara frequencies {ujn = 2m:T) of the func- 
tion fiiUJn) 



r'E/(*^") 



1 

27ri 



/(^)- 



1 



,(iz _Y 



dz. 



where the contour Cmats is given in Fig. [TOl If 

f{z)NB{z) —J- when \z\ -^ oo then it is possible to close 
the contour by adding the arcs of C„„,„; „■ „. in which 

■^ *— ' SGiiii — circ ' 



case C„ 



r<+ 

semi— circ 



a 



semi— circ' 



In general /(z) will 



only have poles or branch cuts on the real axis, such that 
the contours can be contracted along the real axis, which 
gives 

1 r°° 

-- f [1 + NB{uj)Mf{-uj + iQ)]diJ, 

TT Jo 

(Dl) 



where 2i3[/(w + iO)] = /(w + iO) - /(cj - iO) and it 
was used that w < for the contour C~^^^_^j.^, such 
that it is more convenient to take w — > — cj and use 
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Im 



C" 



z] 



^--mats 
A ^ semi-circ 



Re[z] 



FIG. 10: The contour for the bosonic Matsubara summation, 
where Cmats along the imaginary axis is used to rewrite the 
frequency summation as a contour integral. The contours 
CsJmi-circ ^^^ ^^^ posltive and negative energy semicircular 
contours obtained by adding arcs at infinity to Cmats- 



Im[z] 

A 

^ semi-circ ; i ^box 


r 

^mats 
J, ^ semi-circ ; Re[z] 


t 


r 









(^ 



FIG. 11: The same as in Fig. llUl but for the fermionic case, 
where Cbox, the contour between < lj < /i, is related to 
processes inside the Fermi sphere. 



NBi^iLi) = — 1 — Nsiio) ■ The first line can be inter- 
preted as due to thermal gluons and the second line due 
to thermal (anti)gluons and a vacuum contribution. The 
fact that the gluon is its own antiparticle will be reflected 
by3[/(-u; + zO)] = -5[/(^ + iO)]. 

Let us evaluate the above for the specific case f{z) ~ 
\ii[g{z)] which satisfies 5[(7(— a; + iO)] = —'^[g{oJ + iO)]. 
Expanding the logarithm in terms of its real and imagi- 
nary parts 



ln[g(z)] = In \g{z)\ + i arctan 



/S[5(^)] 



nre{-di[g{z)])sgn{^[g{z)]) 



then Eq. (|D1[) can be written as 
/3-i^ln[g(*c.„)] 
1 



duj{l + 2NB{uj)) 



arctan 



S[g(^+)] 



~7re(-5i[.9(^+)])sgn(3[.9(^+)] 



which will be used in the calculation of the RPA correc- 
tion to the thermodynamic potential. 

For fermionic Matsubara frequencies a similar deriva- 
tion can be done, but now the chemical potential is in- 
cluded by writing /(iw„ + /i) 



r'E^(*^- 



-/i) 



1 



/(^) 



1 



3/3(z-m) + 1 



dz, 



where the contour Cmats is given in Fig. [TT] and N{z) = 
(exp(/3z) -I- 1)~^ has poles at the fermionic Matsubara 



frequencies with residue —1. li f{z)N{z) — > when \z\ — >■ 
(X), then it is possible to close the contour as given in 
Fig. [11] by adding the arcs of Cg^^^_^i^^ and the lower- 
and upper-boundaries of Cbox- The full Matsubara sum 
can thus be written as a contour integral over Cmats = 

Cstmi-circ + C'box + C's^mi-circ' ^gaiu if f{z) Only haS 

poles and branch cuts on the real axis, contracting the 
contours along the real axis gives 



N{lo - M)S[/(a; + iO)]da; 
[1 - N{fi - u})]^[f{uj + iO)]duj 
[l-N{uj + n)]^f{-~uj + iQ)]duj. (D2) 



1 ri^ 



TT Jo 



where in the first-to-last line it was used that < w < /i 
for the contour Cbox, such that it is more convenient to 
write N{uj - ^) = 1 - N{fi - w). Similarly for C~^^_^-^^.^, 
where w < 0, we take uj — > — cj in the integral and write 
N{—uj—^) = l — N{u!+iJ,). From the above it is clear that 
the first line corresponds to particles above the Fermi 
sphere, the second line to those in the Fermi sphere and 
the last line is due to anti-particles and contains a vac- 
uum contribution. Note that the above formula must 
be used with care, since one has to check if the above 
simplifications are valid on a case-to-case basis. 
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Appendix E: The gluon self-energy 



where Ff and Gf can be obtained by 



The gluon self-energy, or the so-called polarization ten- 
sor, is given to the lowest order by 



yrab 



9'^A I ^Tr,h,Gof{P + Qh.GofiP)] 

(El) 



(5"''n 



f.i^'^- 



The polarization tensor is purely transverse, i.e., 
QfJ^f,fj.i^ = 0, which can be easily seen using 

i(^^G,hP + Q)-G-hP), 



and the cyclicity of the trace. Expanding the polarization 
tensor in terms of the longitudinal and transverse parts 
relative to the three-momentum q gives 



n/^p.^^/^^'^ + G/P, 



flV "> 



Ff = 



Q' 



n 



/.oo, 



G/ = X (n/,MM - Ff) ■ 



(E2) 



Noticing that even in the imbalanced case, the contribu- 
tion from each flavor simply adds up to the total polar- 
ization tensor, therefore, throughout this section, we will 
not explicitly write out the subscript /, and all the formu- 
lae presented here are applicable to any flavor. However, 
in the main text, we use 11^,^, F and G as the summation 
over all flavors if the subscript / is not explicitly shown. 

After performing the Matsubara sum, the individual 
processes involved can be identified 






' 2 



(27r)3 



E 

Si,S2=±l 



iWg-eo(p + q) + eo(p) 
J2ss' [^~s iPht^Us' (p + q)us' (p + q)7^Ws (p)] 



[N{p + q)(l - N{p)) - (1 - 7V(p + q))7V(p)] 
[iV(p + q)iV(p) - (1 - N{p + q))(l - iV(p))] 



iujq - eoip + q) - eo(p) 

E...[^..(p)7.^.(p + q)MP + q)7.^.(p)] [ _ ^(^ ^ ^^^^, _ ^^ _ ^^^ ^ ^^ )j 

lUJq + eo[p + q) + eo(p) 
J2ss' [^s (Ph^^vs' (p + q)w<,' (P + qht^Vs (p)] 



iujg + ea^p + q) - eo(p) 
d^p Tr[P^^ (p)7o7p^l^ (P + q)707.] 
(27r)3 iujq - S2eo{p + q) + sieoip) 



4(1 - N{p + q))N{p) - N{p + q)(l - N{p))] 



[N^%p + q){l^N-^^{p))-(l-N^'^{p + q))N^^{p)] 



(E3) 



where iV*(p) 



cxp[/J(seo(p)-/j) + l] 

the Fermi distribution of the corresponding flavor and 



iV+(p) = N{p) is where the various functions have been defined as 



Nip) 



cxp[/3(£o(p)+A')] + l 



1 — A^ (p) is the Fermi dis- 



tribution for the corresponding anti-particle. The first 
line is related to particle-hole creation and annihilation, 
the second and third line to particle-antiparticle (and 
holc-antihole) production and annihilation, and the last 
line to antiparticle-antiparticlehole creation and annihi- 
lation. Note that, as expected, the first line vanishes in 
the T, /x — >■ limit and the last line vanishes in the T — > 
limit. The second and third line contain infinite vacuum 
(T, n — 0) contributions, which need to be renormalized. 
The gluon self-energy will be decomposed in terms of a 
matter part, a renormalized vacuum and an infinite vac- 
uum contribution according to 



Tl,,,{Q,T,^l) 



TTinat 



Trrcn.vac. 



TTinf.vac. 



Tjmat 



n^,{Q,T,^l)-^Ilxx{Q,o,o), 



w 



p 

rcn.vac. M^ 



ninf.vac. ^ P^Q 



nAA(Q,o,o)-Q2 
nAA(Q,o,o) 



nAA(g,o,o) 



Q2^0 



flU 



Q^ 



For the vacuum expressions it is possible to extract the 
projection matrix P^w, since the polarization tensor is 
purely transverse and for T, fi ~ there is no preferen- 
tial frame such that the tensor is built up out of only two 
possible quantities (5^,^ and Q^Qu- To include the dy- 
namical properties of the vacuum expression a renormal- 
ization procedure is necessary. In the above the renor- 
malized vacuum expression was obtained by extracting 
the infinite contributions from the vacuum polarization 
tensor at the renomalization point Q^ = 0. The reason 
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for extracting the factor Q^ from the vacuum IIaa before 
setting (5^ = is because in the vacuum the dressed glu- 
ons are massless, i.e., IIaa oc Q^, c.f. Eq. (jE7p . Therefore 
Q^ = is still a pole for the vacuum expression, at which 
the residue is 1 after the renormalization. In conclusion, 
the following full renormalized polarization tensor will be 
used 



In the high-density limit (/x ^ m), a good approxima- 
tion is to evaluate the matter part for the case of massless 
quarks. Thus for massless quarks in the zero-temperature 
limit the matter parts of IIoo and 11^^ can be calculated 
explicitly 



jjren 



n„ 



n' 



inf. vac. 



Trmat 



Trren.vac. 



(E4) 



n™^*(r = o) 



27r2 



n™;*(T = o) 



1 

24 



E 

Si,S2=±l 



lS\LOq 



(Q^ + 2(j2 



+ E 




1 + 



S2 = ±l 



q+iUq 



In 



1 



2siS2/i 



S2 



q - ISlUJq 

H{3Q^ - 4^^) ^^ iujq 



r 



24g 

Si ,S2— ±1 



4 q 
iSlLOq) In 



In 



1 



2S2/^ ■ 
q—iuq 



2s2t^ 



q+iujg . 
q - iSlLOq_ 



(E5) 
(E6) 



These expressions give the same result as can be found 
in Refs. [40l - |42| . Performing an analytic continuation 
to real time {iujq — )■ w -f iO), the above expressions have 
branch cuts from max(0,g — 2/i) < |ci;| < g and q < 
\uj\ < q + 2/i, the origins of which can be found from 
Eq. ()E3P and are due to particle-hole processes and finite- 
T, /i contributions to particle-antiparticle production and 
annihilation, respectively. Because Lorentz invariance is 
broken due to the presence of the gas, n™'** is a function 
of go = i{i<jJq) and q separately since it can be a function 
of qfj,Ufj, and q = a/Q^ - (g^u^)^, where u^ = (1,0,0,0) 
defines the rest frame of the system [l^ . 

The renormalized vacuum part 11'^™™'^' has been found 
inRef. \M 



P O'^ — 



d.Ta;fl — x) In 







1 



3^"^ 4^2 



- 1- 



m^ + x{l — x)Q'^ 
2to2 



X 1 



1 + ^^ArcCothJl 



Q^ 



Q^ 



2n2 



^P - 



g'Q 



f+ln 



i^) 



247r2 



(E7) 



where it was assumed that the quarks have equal masses 
and in the last line it was expanded for small masses. 
In real time this expression has a branch cut for |a;| > 
■y/g^ + 47r? ~ q, which is due to particle-antiparticle pro- 
duction and annihilation. 

Using a Sommerfeld expansion temperature correc- 
tions to the gluon self-energy can be obtained. The Som- 
merfeld expansion can be summarized as 



Niep - At) =&i^i - ep) - E ^(1 - ^'-'nCi^n) 

a2n-l 



X 



d^'^'^ 



m)t 



■2n 



(E8) 



The derivation is analo gou s to that of the non-relativistic 
Sommerfeld expansion |44l | . One has to keep in mind that 
the above is not a complete expansion for small temper- 
atures since also the chemical potential depends on tem- 
perature. In this manner the T^ correction to the gluon 
self-energy is found to be 



n 



00 



TTinat 
^^00 



(T: 



0) =^ir2 y . 

2 6 ^^ 

s=±l 



2s/i 



n" 



9' 



1, 



TTinat /Ti r\\ y _l7^2 

t-t- ^ -''>- 2 6 [{ioj, 



2q 



In 



(iujq — q){ii^q + 2s/i + q) 



{iujq + q){iijq 
-(2/i-g)(2^ 



' 2s/i 



q) 



(2/.-g)2][(zc.,)2-(2M + g)2 



(E9) 



20 



Using Eqs. (|E2IIE9|) the behavior of F and G can be 
found in several hmits. Keeping the ratio x — q/uj fixed 
and expanding up to zeroth order in oj, the hard dense 
and hard thermal loop (HDL/HTL) expressions are re- 
obtained [il, [11 
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where 
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is the gluon thermal mass. The reason that the above 
limit returns the HDL and HTL expressions, is due to the 
fact that the HDL/HTL approximation (m = and cj <C 
q <C /^) takes into account only the low-energy processes 
around the Fermi surface, namely particle-hole processes. 
Some other useful limits are 
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where it can be seen that in the static long-wavelength 
limit {u!,q — >■ 0) the electric gluons (F) are screened, 
while magnetic gluons (G) are dynamically screened. 
Furthermore, the large momenta and frequency behav- 
ior of the matter and vacuum parts are separately seen 
to be 
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The vacuum expressions are proportional to Q^ since in 
the vacuum the gluons remain massless and thus Q^ = 
is still a pole of the propagator. The logarithmic Q^ 
dependence is usually absorbed into the vertex and sub- 
sequently interpreted as the varying of the coupling con- 
stant with the energy scale Q. 



Appendix F: The thermodynamic potential of ideal 

gases 

Consider first the well-known ideal Fermi gas term of 
a single species 
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Using Eq. (|D2p with /(z) = ln[— z^ -f- eo(p)] and using 
the principle value logarithm with a branch cut on the 
negative real axis 

lim3[/(a; -f i??)] ^ Iim5[ln(-w2 + e\{^) - 2zcli?/)] 

= 7re[a;2-eo(p)]sgn(-w), 

the ideal Fermi gas contribution becomes 
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In the above it was assumed that the arc at infinity van- 
ishes, which is the case if the time-ordering in the path- 
integral is taken into account properly by multiplying the 
integrand by e"*'^"'' and taking the limit 77 J, 0. The infi- 
nite vacuum contribution is clearly visible and should be 
subtracted, giving for massless fermions 
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The ideal Bose gas follows from 
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which will lead to the same expression in both the Lorentz 
and Coulomb gauges. In the following, the Coulomb 
gauge is taken (/^ = (0, V)^), 
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The unphysical degrees of freedom elearly drop out due 
to the ghost contribution 9^/^. Using Eq. (|Dip with 
f{z) — ln(— z^ + q^), the ideal Bose gas contribution is 



PJ (27r) 
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dui. 



Again the vacuum term is clearly present and will be 
subtracted, giving the Stefan-Boltzmann law 
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Appendix G: The quark self-energy 



For completeness also the quark self-energy will be de- 
rived in the Lorentz gauge to lowest order. 
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Performing the Matsubara sum gives 
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In the last line the indicator function Is2,s3 = Ss2+(^S3+ — 
5s2-Ss3- is exactly the vacuum contribution such 
that AS is the vacuum subtracted self-energy, i.e., 
AE (T = /x = 0) = 0. To first order the self-energy for 
a single quark and antiquark can be found by using the 
free quark energy projectors and helicity projectors 

K'{iojp,p) = Tr[7'^(p)P^^7oS(*Wp,p)]. 

After projecting and evaluating the self-energy at iujp = 
5ieo(p) — 1^ the function in the first square brackets in 



Eq. (|G1|) reduces to J^si,s2,s3{p,<i) defined in Eq. ^ 
In the zero-temperature limit the self-energy for a single 
quark then reduces to 
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Note that the final result is independent of helicity s. 
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